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The exact relativistic treatment of a rotating disk surrounded by a magnetized material halo is 
presented. The features of the halo and disk are described by the distributional energy-momentum 
tensor of a general fluid in canonical form. All the relevant quantities and the metric and electro¬ 
magnetic potentials are exactly determined by an arbitrary harmonic function only. For instance, 
the generalized Kuzmin-disk potential is used. The particular class of solutions obtained is asymp¬ 
totically flat and satisfies all the energy conditions. Moreover, the motion of a charged particle on 
the halo is described. As far as we know, this is the first relativistic model describing analytically 
the magnetized halo of a rotating disk. 


I. INTRODUCTORY REMARKS 

In the observational context, many ambiguities still exist about the main constituents, geometry and dynamics 
(thermodynamics) of the galactic disk-haloes. However, there are several different observations which probe the 
galactic and surrounding galactic magnetic field. A current revision of the status of our knowledge about the magnetic 
fields in our Milky Way and in nearby star-forming galaxies is summarized in [Tj . Additionally, a study of the disk and 
halo rotation are reported in [2], whereas the possibility of magnetic fields can be generated in the outskirts of disks 
is studied in [3]. Solutions for the Einstein and Einstein-Maxwell Field Equations which are consistently applicable 
to the context of astrophysics remains a topical problem. Nevertheless, the effects of magnetic fields on the physical 
processes in galaxies and their disk-halo interaction have been scarcely considered in the past. Similarly, the relevance 
of relativistic models of disks around black holes in a magnetic field is discussed in [4]. 

The presence of the electric field on the dark matter halo models has been considered in [5; , whereas the presence of 
electromagnetic field in the halo-disk system has been studied in mm in which the gravitational sources are statics. 
In [B] we provided an detailed overview of the research in the relativistic disks, accordingly we shall not repeat them 
here. 

In this paper we considered the conventional treatment of galaxies modelled as a stationary thin disk and, corre¬ 
spondingly, we associate the halo with the region surrounding the disk. We present the conformastationary version 
of the static thin disk-halo systems studied in [6]. In addition, we study the features of the principal quantities 
characterizing the dynamic of the magnetized haloes corresponding to the disks presented in [5j. Therefore, we take 
the definition in Ref. [9] as standard, following the original terminology by Synge m- conformastationary are those 
stationary spacetimes with a conformally flat space of orbits. 

Accordingly, we show that the rotating disk-haloes with isotropic pressure, stress tensor and heat flow generalize 
the static disk-haloes obtained in !6j. Our results are compatible with those presented in [5] on possible features of 
galactic halo. Moreover, the description of the motion of charged particles on disk is deduced and is in agreement with 
the results of the similar analysis discussed in m- As far as we know, this is the first relativistic model describing 
analytically the relativistic magnetized halo of a rotating disk. 

The paper is organized as follows. In Section[TTJ the distributional Einstein-Maxwell equations for haloes surrounding 
thin disks are obtained. In Section |Hl| we obtain expressions, in terms of an arbitrary harmonic function, for the most 
important quantities characterizing the dynamic of the disk and halo. In Section |IV| we first calculate quantities for 
an harmonic function described by the generalized Kuzmin-disk potential. Then, we analyze the obtained results and 
calculate the constants of motion of the disk. Moreover, the description of the motion of a charged particle on the 


Finally, we complete the paper with a discussion of the results in Section [VII 


halo is shown in Section Wl whereas the stability of the orbits of test particles in the halo is studied in Section VI 
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II. EXACT SOLUTIONS FOR RELATIVISTIC MAGNETIZED HALOES SURROUNDING THIN DISKS 


In this section we consider the conventional treatment of rotating galaxies modelled as a stationary thin disk and, 
correspondingly, we associate the magnetized halo with the region surrounding the disk. To do so, we formulate the 
distributional Einstein-Maxwell field equations assuming axial symmetry [12]. We also suppose that the derivatives 
of the metric and electromagnetic potential across the disk space-like hyper-surface are discontinuous. To formulate 
the corresponding distributional form of the Einstein-Maxwell field equations, we introduce the usual cylindrical 
coordinates x a = (t, r, z, p>) and assume that there exists an infinitesimally thin disk located at the hypersurface 
z = 0. Accordingly, we identify the halo surrounding the disk with the positive (z > 0) and negative (z < 0) regions 
around the equatorial plane z = 0, denoted here by the superscripts “±”. So that the metric and the electromagnetic 
potential can be written, respectively, as g a p = g+p0(z) + g~p{ 1 — 9(z)} and A a = A+9(z) + A~{ 1 — 9(z)}. Here 9{z) 
denotes the Heaviside distribution. As a consequence, the Ricci tensor reads 

R«p = B+ P 0(z) + R~p{ 1 - 9{z)} + H a pS(z), (1) 


where S(z) is the Dirac distribution and 


R a/3 2 


J z p 




( 2 ) 


with 7 a p = 2 g a p tZ and all the quantities are evaluated at z = 0 + . In agreement with (jTJ) the energy-momentum tensor 
and the electric current density acquire the form 


T a p = T+ 0 9(z) + T~p{l - 0(z)} + Q a p5(z), 


(3a) 


J a = Ja9(z) + J a i l ~ 9{z)} +I a S(z), 


(3b) 


where and are, respectively, the energy-momentum tensors and electric current density of halo. Moreover, 
Q a p and X a represent the par t of the energy-momentum tensor and the electric current density of disk. The energy- 
momentum tensor in (3a) is taken to be the sum of two distributional components, the purely electromagnetic 
(trace-free) part and a “material” (trace) part, 




( 4 ) 


where E±p is the electromagnetic energy-momentum tensor 


= FavFp v — -g a pF liI/ F lil ', (5) 

with F a p = Ap tC[ — A a fi and M±p is an unknown “material” energy-momentum tensor to be obtained. Accordingly, 
the Einstein-Maxwell equations, in geometrized units such that c = 8irG = go = eo = 1, are equivalent to the system 
of equations 


R % - \g<*t3 R± = F% 3 + M±3 

( 6 a) 

R a(3 ~^9a/3 R — Qck/3? 

( 6 b) 

TZ a 0 _ T a 

b ± ;/3 — J ±i 

( 6 c) 

[F a P] n p = I a , 

( 6 d) 


where H = g a @H a p. The square brackets in expressions such as [ F a/9 ] denote the jump of F a P across of the surface 
z = 0 and n p denotes a unitary vector in the direction normal to it. To obtain a solution of the distributional 
Einstein-Maxwell describing a system composed by a magnetized halo surrounding a rotating thin disk we shall 
restrict ourselves to the case where the electric potential A t = 0. 

To solve the Einstein-Maxwell equation § we assume the conformastationary space-time given by the line element 

ds 2 = -e 2cl> (dt + todip) 2 + e~ 2 ^{dr 2 + dz 2 + r 2 dp 2 ), (7) 


where <f> depending only on r and z and fd is an arbitrary real constant. We also assume that the magnetic potential 
A v is time-independent. Accordingly, by computing the Einstein tensor G a p explicitly from the line element (J7| and 
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the electromagnetic energy-momentum tensor ([5]) , we obtain for the non-zero components of the energy-momentum 
tensor of the halo (EMTH) M±p = E±p - : 


= - e 2 ( 1+/3 ^{/? 2 V0 • V(t> - 2/3V 2 <£ + ^- 2 e 2W VA v • VA v - | r - 2 e 2 ( 1+/, )*Vw ■ Vw} 

= e 2(1 +W{ |Vw ■ V(t> + |jr- 2 e 2(1+W wVw • Vu - /3 2 wV</> • V</> + 2/3wV4 + ^Vw • V0 

- |r 2 e 2W wVi», • + ^V 2 w - r _1 Vw ■ Vr} 

M± = (1 - /3)V 2 0 - (1 - /3)<4, + (/3 2 - 2/3)4 + 4> 2 Z - \r~ 2 e 2 ^{Al tr - A 2 , z ) 

+ l r -2 e 2(l+W (w 2 _ w 2) 

M tz = ^ _2 e 2(1+W w r w ~ - (1 - /3 2 + 2/3) <4 4 - (1 - /3)<4 Z - r^e 2 ^ A^^A^^, 

Mf z = ~r- 2 e 2 ( 1+ W(u; 2 r - u, 2 ) + </> 2 r - (1 - /3)0, zz + (1 - /3)V 2 <(> + (/3 2 - 2/3)4 

I I,,-2p 2 / 3 0('/l 2 _ /l 2 / 

' 2' e 

= r 2 V<(> • V0 + (1 - /3)r 2 X7 2 (f> - (1 - /3)rV</> • Vr - Je 2W V^ • VA^ 

+ e 2(1+w { I(i + 3 r- 2 e 2W w 2 )Vw • Vw - /3 2 w 2 V</> • V</> + 2/3w 2 V 2 <)> + wV 2 w 

- 2r -1 wVw • Vr + (3 + /3)uVw • V<£ - ^"V^^VA^ • VA^j. 


Moreover, from (6c I the non-zero components of the electric current density on the halo have the form 

4 = e - (1 - 3W V • {ur" 2 e (1+W V4}, 

JV = e -(l-3/3)0 v . | r -2 e (l+/3)^ v ^| ] 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 


(8f) 


(9a) 

(9b) 


where all the quantities depending on r and 2. 

The discontinuity in the 2-direction of Q a p and 1“ defines, respectively, the surface energy-momentum tensor 
(SEMT) and the surface electric current density (SECD) of the disk S a p, more precisely 


*5a/3 — J Qapfi^z'jdSn — yjQzzQafi ; 
J a = Jl a S(z)ds n = ^l a , 


(10a) 

(10b) 


where ds n = yjg zz dz is the “physical measure” of length in the direction normal to the 2 = 0 surface. Accordingly, 


for the metric Q, the non-zero components of S a p and J a are given by 

S tt =4/3e( 2+ «4> z , (Ha) 

S tv = e( 2+ W(4/3o4 z + w,*), (lib) 

S rr = 2(1 - P)e~ p U, z , (He) 

S vv = e ( 2 +/ 3 )'3{(4 poo 2 + 2(1 - 4^ 2 e _2(1+/3)0 )0,z + 2 ww z }, (lid) 


t _ _-2 . 3/90 


T =r 


c r ~' r u) [A^, jZ ], 


- „- 2 p3/3</> 


J* = -r 


[4, 




(12a) 

(12b) 


respectively. Note that in (111 and (12) all the quantities are evaluated on the surface of the disk (2 = 0). 

In order to reduce the complexity of the last field equation systems we assume that the halo’s electric current 
density vanishes (i.e. J± = J± = 0 in (|9|)), it turns out that the magnetic potential and the metric functions (j) and 
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u! become completely determined in terms of an arbitrary harmonic function U ( r , z) as follows (see [B, ;Bj for more 
details), 




- \ TV 

(13a) 

l ru - 

(13b) 

— In (1 — U), 

(13c) 

fouj U j 

(13d) 


with k and k u arbitrary constants. Since the non-zero components of the EMTD and EMTH and the electric 
current density directly depend on the metric functions and magnetic potential, we observe that the entire solution 
is determined by a single harmonic function. 


III. EXACT RELATIVISTIC MODEL FOR MAGNETIZED DISK-HALOES 

So far, by using the inverse method and the distributional formulation of the Einstein-Maxwell equations, we have 
obtained the separate energy-momentum tensors of the disk and halo. In addition, we have discussed out a method 
to determinate it in terms of an arbitrary harmonic function. Now, the behavior of the energy-momentum tensors 
obtained must be investigated to find what conditions must be imposed on the solutions and the parameters that 
appear in the disk-haloes models in such a way that it can describe reasonably physical sources. We shall now study 
the possible features of the disk by assuming that it is possible to express its energy-momentum tensor in the canonical 
form 


S a p — (m + P)V a V /3 + Pg a /3 + QaVp + QpV a + n Q/3 , (14) 

where Q a V a = Q a V a = 0, a = ( t,r,ip) and all the quantities are evaluated in z = 0 + . Similarly, we assume that its 
possible to express the energy-momentum tensor of the halo in the canonical form 

= (^ + P^VaVp + p ± g a p + Q±Vp + Q$V a + n± (15) 


where Q ± a V a = Q ±a V a = 0, a = ( t , r, z, <p) and all the quantities depend on r and z. Consequently, we can say that 
the disk and halo are constituted by a some mass-energy distributions described by the energy-momentum tensors 
(14) and (15), respectively. The V a is the four velocity of certain observer. Correspondingly, /z, P, Q a and n a/ 3 
are then the energy density, the isotropic pressure, the heat flux and the anisotropic tensor on the surface of the 
disk. Analogously, /i ± , P ± , Q± and rp, 3 are then the energy density, the isotropic pressure, the heat flux and the 
anisotropic tensor on the halo, respectively. Thus, it is straightforward to see that for the halo we have 


= M± 0 V a VP, (16a) 

P ± = \u a ?M%, (16b) 

2a = ~^V a - M%V^ (16c) 

= HSHfKM* - P ± P M ,), (16d) 


where the projection tensor is defined by TL^ = g^ + V), V v and all the quantities depending on r and z. Whereas , 
for the disk we have 


n = S a pV a VP, (17a) 

P = \u ap S a p, (17b) 

Qa = -fJ.V a - S a pV (17c) 
n«/3 = - PH„ V ), (17d) 


where all the quantities are evaluated in z = 0 + . It is easy to note that by choosing the angular velocity to be zero in 
Equation (A7 ) we have then a fluid comoving in our coordinates system. Hence, we may introduce a suitable reference 
frame in terms of the local observers tetrad (A3) and (|A4|) in the form 




(18) 
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with the corresponding dual tetrad 


{V a ,I a ,K a ,Y a } = {-hW a ,h< r l,hW,hM}. 


(19) 


Since the SECD of the disk J a can be also written in the canonical form J a = aV a + jY a , a can be interpreted 
as the surface electric charge density and j as the “current of magnetization” of the disk. A direct calculation 
shows that the surface electric charge density a = 0, whereas the “current of magnetization” of the disk is given by 
j = — r -1 e 2 ^ where, as above, [A VtZ \ denotes the jump of the z—derivative of the magnetic potential across 

of the disk and, all quantities are evaluated on the disk. 

By using the results obtained in the precedent section, we can write the surface energy density of the disk and the 
energy density of the halo can written as 


MO) 


4/317, z 

(/3 + l)(l-f/)w 1 


( 20 ) 


and 


,*<,.)={<*+- m ^ 1 } . <*> 

respectively. Moreover, we have a barotropic equation of state on the surface of the disk, which can be given by 
P(r ) = iyj, , with 77 = (1 — /?)/3/3, in such a way that we can use the energy conditions and the causality requirement 
for the speed of sound on the disk to obtain the physical range of possible values of the parameter /3. Analogously, 
the pressure of the halo we have P ± (r, z) = (r, z ), where 

(4 — 2/3 — /3 2 )r 2 — r2e ~ 2< ^ + k ^^' > (l + 3fc 2 r _2 t/ 2 e 2/3<i4 (l — e 2| 3)) 

3 ((2/3 + ^y-^rV 2 H^Mj 

in such a way that the pressure of the halo not only depends on the energy density but also on the gravitational and 
magnetic fields through the function 0. The Heat function of the disk is given by 


Qa(r) = 


( 22 ) 


Similarly, by inserting (15) into (16) we obtain for the heat flux of the halo 

Le( 1+2 «« 


Q = 


2(1 + /3)r 


{2(1 + /3 )U r - (3 + f3)r(U 2 r + [7 2 ) e ( 1+/3 ^} 6%. 


(23) 


The non-zero components of the anisotropic tensor of the disk read n vv (r) = r 2 H rr (r) where 


2(1 - /3 )U Z 


n rr (r) — ^ 

3(1+/3)(1 -U)—e 

Moreover, it is easy to see that the anisotropic tensor of the halo reads 


(24) 


T± 

a/3 


o±T j 

r cc. 1 


P±K a K p + P±Y a Y 0 + 2P±I [a K 0) 

(25) 

p± _ p 2/30TT± 

J- ly* O' A.A.j-.ty, * 

(26) 

Pt = e 2 ^n±, 

(27) 

, e 2 ^ 


p ± = n* 

¥> r 2 VV' 

(28) 

P± = e 2/3 ^n± . 

(29) 


where 
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and 


_ e 2(i+^ f / k 2 (l + P ) 2 2(l+/3) : ~ _ 2 ^ _ 

J-J-rr r* /-. . /TV O O I I "T” 1 O 70 


3(1 +/?) 2 r 2 l V 2 
+ (-/£(! +/3) 2 - 


fc 2 


4(1 + p - p 2 )r 2 - 3 ^( 1 + ^) 2 . r - 2 [/2 e 2W (1 _ e 20^ V 2 


2 1 1 , o\2 (1 + /^) 2 2 - 


20 


fc 2 


+ 2(1 + p - p 2 )r 2 - 3fc - (1 4 +^ 2 r- 2 t/2 e 2^ (1 _ e 2 <^ u* 


-3(l-0 2 )rV< 1+ ^l7 )rr j, 


(30a) 


e 2(l+/3)0 

n ± =- 

22 3(1 + P) 2 r 2 


-ki{i + py- 


2 (1+/3) 2 9 _ 


/c 2 


-r e 


20 


n w = 


+ 2(1 + P - P 2 )r 2 - 3 ^ (1 4 + ^ )2 r-2[/2 e 2W( 1 - e 2 ^)) U 2 
+ + ?il±^l! r 2 e -20 _ 4(1 + p _ p2y _ Zkt(l+P) 2 r - 2£/ 2 e 2/30 (1 _ g 20^ p-2 

-3(l-/3 2 )r 2 e-( 1+ Wf/, 2z |, 

^^|^^{2(1 + P - /3 2 )r 2 - ii±^! r 2 e- 2 ^ + (1 + 3 A£r" W**(l - e 2 *)) } 


(30b) 


(1-/3) 
1 + P 


re (1+w /7, 


e 2(l+/30) 

n± = - 

r2 (1 + / 3) 2 


| ( - 2(1 + p - p 2 ) + (1 + J )2 + fc “ (1 2 ^ )2 )- (1 - /3 2 )e-( 1+ ^[/ rz J. 


(30c) 


(30d) 


Notice that 'P ± = P* + Pp + P± = 0 and, consequently, the trace n ±a a = 0. We have obtained expressions for the 
energy, pressure and the other quantities characterizing the dynamic of the halo. All the dynamic quantities have 
been expressed in terms of an arbitrary U(r,z) harmonic function. Finally, as we know, the electric current density 
of the halo is zero whereas it is easy to note that the magnetization current density on surface of the disk is 

[Ur] 


j(r) = - 


, s 3/3 

k( 1 - C/)!+3 


(31) 


It is important to remark that k u is a defining constant in (221 and (231. Indeed, when k w = 0 the heat flux functions 
Q a and Q* vanish, a feature of the static systems. Due to we used the inverse method, no “a priori” restriction are 
imposed on the physical properties of the material constituting the disk and halo. The non-zero components of the 
energy-momentum tensors of the disk and halo result of “the nature” of the chosen metric and the corresponding 
solutions. So, in our case, the non-zero component S rr and Sup are conditioned by the parameter p and the metric 
function u in such a way that when p = 1 the component S rr vanishes, whereas Sun = 0 when ui vanishes. The 
decomposition of the energy-momentum tensor of the disk-halo system into (14) and (15) were chosen with the aim to 


describe the SEMT and EMTH by the more general fluid model. Hence, the heat flux appear here in a “natural” way 
as a function determined by the metric function u and, consequently, by the “rotation”. Unfortunately, as we can see 
from (17c) and (16c), this function is oriented along the closed circular orbits and thus its physical interpretation is 


unclear. It is an issue that remains unanswered in this manuscript, but should be addressed in the future. 


IV. ROTATING KUZMIN-LIKE DISK WITH MAGNETIZED HALOES 


As an example of application of the formalism described in the precedent sections, we now consider the magnetized 
haloes surrounding the rotating disks generated by a generalization of the Kuzmin-disk potential in the form mm 

Nwz/r) *«/*) = <- 


U = -J2 


n —0 


Rn+l 


1! dz n 


(32) 


where P„ = P n (z/R) is the Legendre polynomials in cylindrical coordinates that was derived in the present form by 
a direct comparison of the Legendre polynomial expansion of the generating function with a Taylor of 1 /R 1141. b eing 
R 2 = r 2 + z 2 and b n arbitrary constant coefficients. The corresponding magnetic potential, obtained from (h 3 ab and 


(13b), is 


N 


A 


= Jw 5 (zir 

k 2^ ° U 77.1 


n—0 


dz 7 


d) 


(33) 
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(a) ( b ) 

FIG. 1. Surface plots of the energy density (a) and (b) p* on the halo as a functions depending on r and z with parameters 
a = bo = &i = k = k u = 1 and /3 = 0.75. 


where, we have imposed +1^,(0, z) = 0 in order to preserve the regularity of the axis of symmetry. Next, to introduce the 
corresponding discontinuity in the first-order derivatives of the metric potential and the magnetic potential required 
to define the disk we perform the transformation z —> \z\ + a. It is worth noting that, for the two first members of 
the family of solutions (N = 0 and N = 1) we have 


and 


U 0 = - 



k 

s/f 2 + {\z\ + l) 2 ’ 

fr 0 (|z| + 1) 

k^Jf 2 + (\z\ + l) 2 ’ 


(34a) 

(34b) 


-Api — — 


U! = ~ 


k(\z\ + !) 


1 - 


ksjr 2 + {\z\ + l) 2 1 & 0 (|2| + 1) ((|z| + l) 2 + f 2 ) 


1 + 




V'r 2 + (|z| + l) 2 I 6 0 (( \~z\ + l) 2 + f 2 )\' 


(35a) 

(35b) 


respectively, where b 0 = b 0 /a and b 1 — b^/o? whereas f = r/a and z = z/a, moreover A v = A v /a. For the two first 
members of the family of solutions the surface energy density of the disks can be written as 


and 


Mo 


_ ^k _ 

/_ ._ x (2/3+l)/(/3+l) ’ 

(1 + p){r 2 + 1)03+2)7(2/3+2) ( bo + TFTl) 


Mi = 


4/3{(6 0 -k)r 2 + k + 2 M 


(1 + /3){r 2 + l)(2-/3)7(2/3+2) j ( f 2 + 1)3/2 + ^ 2 + X ) + ^ } 


(2/3+l)7(/3+l) • 


(36) 


(37) 



















respectively. 


Similarly, for the two first members of the family we have the heat flux of the disks 


Q aO 


Q al 


b 0 k ^a 

Vf 2 + 1 + b 0 ’ 

(b 0 (f 2 + l) + 6.) 

(f 2 + l) 3 / 2 + b 0 (r 2 + 1) + b 1 ' 


and the corresponding anisotropic tensor 


_ 2 ( 1 -/ 3 ) 6 0 _ } 

3(1 + /?) + & 0 ) 1/(1+/3) (f 2 + l)(3/3+ 2 >/(2+ 2 ^) 

2(1-/?) ((60-^ + 60+26.) 

3(1 + /3) {(f 2 + l) 3 / 2 + b Q (r 2 + 1) + b^ 1/(1+f>) + l) (5/3+2)/(2+2/3) 


(38a) 

(38b) 


(39a) 


(39b) 


As we know, the another quantities are P = (1 — /3)/z/(3/3) and II™ = r 2 H rr . In the last expressions we have used the 
dimensionless expressions jl = a/i, II W = all rr Finally, for the two first members of the family we have the current 
of magnetization as 


/_ x 20 /(i 3 +i) ’ (40a) 

k(f 2 + l)(3+/3)/ ( 2 + 2 / 3 ) U + v^TlJ 

2r (b 0 {r 2 + 1) + 3^) 

— r I - x 20/(i+0 )' (40b) 

k(f 2 + l)(5-0)/(2+20) I ( f 2 + 1 ) 3/2 + + 1) +& 1 J 
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where j = aj and we first have assumed that the ^-derivative of the magnetic potential present a finite discontinuity 
through the disk. In fact, as we have said above, the derivatives of U and A v are continuous functions across of the 
surface of the disk. We artificially introduce the discontinuity through the transformation z —> \z\ + a . 

It is worth noticing that the mass surface density as well as the isotropic pressure of the disk decay very rapidly 
(as 1/r 3 and I/ 7' 5 for N = 0 and N = 1, respectively) indicating that the above solution can be associated with a 
disk with a finite energy-momentum distribution. In every case, the characteristic size can be adjusted through the 
parameters and 61 of the solutions. Moreover, a simple calculation of the curvature invariants reveals that the 
solution is asymptotically flat and singularity-free. 

To illustrate the results corresponding to the principal quantities describing the halo in Fig. [T] we show the behavior 
of energy densities gA on the halo as a function of r and 2. In each case, we plot /ij(r, z) [Fig. [lja)] and /zf(r, z) 
[Fig. [TJb)] for the indicate values of the parameters. It can be seen that the energy density is everywhere positive 
and vanishes sufficiently fast as r increases. 

In Fig. pi we show the behavior of pressure P ± on the halo as a function of r and z. In each case, we plot Pq (r, z) 
[Fig. j2|(a)] and P±(r,z) [Fig. [2j[b)] for the indicate values of the parameters. We can see that pressure is always 
positive and behaves as the energy density of the halo. Thus, we can see that the behavior of these quantities are 
in agreement with the results published in [Sj. Moreover, we also computed these functions for other values of the 
parameters within the allowed range and in all cases we have found a similar behavior. 


A. The constants of motion 


To proceed further, we evaluate the constants of motion. Therefore, from (|13c|) we have 

1 


In 


1 


1 + 0 \l-UJ ' 


(41) 


Then, for the solution (34a I we may write 


<(>o — 


1 + 1 


In 


R 2 + 2\z\ + 1 


R 2 + 2151 + 1 + b n 


(42) 


where R 2 = r 2 + z 2 . This follows that the metric potentials gtt and (jt v for R —> 00 in the disk (z = 0 ) become 


g tt ~ -1 + 


2 b n 


* 2(3 + 0 , 0 (±_ 


(1 + P)R (1 + P) 2 R 2 \R 3 J ’ 

1 


_ k ul b 0 2kyb 2 
9t<P ~ R (1 + P)R 2+ \R 3 

This implies, as is well known (See [T5]), that the total mass-energy of space-time associated with the disk is 

b„ 


(1 + 0 ) 

On the other hand, in ( x,y,z ) coordinates we find that 

k “b 0 y . n 

QqI — -~- O , 

R 3 VI? 4 


9 02 — 


k u b 0 x ( 1 


R 3 


R 4 


903 — O 


R 4 


(43a) 

(43b) 

(44) 

(45a) 

(45b) 

(45c) 


As an application, we use the same procedure as in m and see that the angular momentum Lmo is in the 2 -direction 
and is given by 

Lmo = ~jk u b 0 . 


( 46 ) 
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According to (34b) the magnetic field is 


Bn = ~ 


b n ? 


(?e r + ^5+ 


(47) 

k(R 2 + 2|5| + 1 ) 3 / 2 ' ' - ' ' 

where e a are unit basis vectors in cylindrical coordinates. Accordingly, by expressing the components of the magnetic 
field in Cartesian coordinates and taking the limit as R —► 00 of Bo(x, y , z ) and by using the formula (44.4) of Landau 
and Lifshitz m we may conclude that the magnetic momentum may be written as 

Lbo = ~r- (48) 

k 

We thus see that constants k and k u defines the gyromagnetic ratio L M0 /L B0 = (kk UI )/2. 


V. MOTION OF A CHARGED TEST PARTICLE IN THE HALO 



FIG. 3. Surface plot of the velocity (a) v 2 and 0 -slices of the surface plot of the velocity (b) on the halo as a functions 
depending on r and 0 with parameters a = bo = 61 = k = = 1 and /3 = 0.75. 


The motion of a test particle of charge e and mass m moving in the halo is described by the following Lagrangian 
density: 


£ = T.ga^V 13 + — A 0 


(49) 


where g af 3 and A a are, respectively, the components of the metric and electromagnetic potential, given here by ( |13| ) . 
The equations of motion of the test particle can be derived from (|49| by using the Euler-Lagrange equation. Then, 

(50) 


The velocity of the particle as measured by the local observers is given by v a = u*(f a + fl</? a ), where 

(1 - [ 7 ) 1 /( 1 +« 


v* = 


(1 + k u Ufl)V^- — v 2 


(51) 
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FIG. 4. Surface plot of the velocity (a) v\ and 2 -slices of the surface plot of the velocity (b) on the halo as a functions depending 
on r and 2 with parameters a = bo = b\ = k = k u = 1 and /3 = 0.75. 


Here, the 3-velocity v and the angular velocity f 1 of the particle as measured by the local observers are given by 


rH(l - U) 

1 + k„un 


(52) 


and 


ME# + U D ((1 + P) + T^u) ± y/(U?r + U%)D 
n = 2(1 + P)r{ 1 - U) 2 U t r - 2(E# + U%)A ’ 

D = 4(1 + p)r{l - U)U r + (E# + U%) (k 2 Jl + /3) 2 - 4/3r 2 ) , 

A = p r 2 (l-U) + klU^l + p+^^j , 


respectively. All the quantities depend on r and 2 . In Fig. |3](a) and Fig. 3](a) we show the behavior of the velocity v 2 o 
and v 2 i of a charged particle following an “magnetogeodesic” motion on the halo for the values of indicated parameters, 
respectively. Additionally, in Fig. gb) and Fig. gb), we plot the 2 -slices of the surface plot of the velocity and 
v 2 q and v 2 \ for the indicated values of the parameters, respectively. These curves are obtained via vertical slices of 
the surface v 2 = u 2 (r, 2 ) (a vertical slice is a curve formed by the intersection of the surface v 2 = v 2 (r,z) with the 
vertical planes). For each curve, we can see that the velocity is always less than 1, its maximum occurs around r = 0, 
and it vanishes sufficiently fast as r increases. It can also be observed that the maximum of the velocity decreases 
as the values of 2 increases. We also computed these functions for other values of the parameters within the allowed 
range and in all cases we found a similar behaviour. Naturally, the description of the motion of charged particles on 
disk here deduced is in agreement with the results of analysis of the magnetogeodesic motion of the particle in the 
magnetized disks discussed in El- 
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FIG. 5. Plot for the variation of V"(r) = for the two first members of the family of solutions (a) Uo and (b) U\ for values 
of the parameters given by /3 = 0.25, L = a = = bo = bi = 1 and E = 50. 


VI. STABILITY OF ORBITS OF PARTICLES IN THE HALO 


Since the Lagrangian density (491 does not depend explicitly on variables t and tp, the following two conserved 
quantities exist 


dC 

Pt = Ti = 

dC 


E 

1 

m 

L 


Plfi — ' — ! 

oip m 


(54a) 

(54b) 


where L and E are, respectively, the angular momentum and energy of the particle as measured by an observer at 
rest at infinity. In the halo near the surface of the disk (z —> 0 + ) the motion equations can be reduced to the form 
f 2 + V(r) = E 2 /m 2 , which describes the motion inside an effective potential V given by 


V(r) = 


43 


m 


V 2 (1-C/)i+3 


E 2 r 2 [(l-U)^ + ^J(1-U)w -(1-U) 2 ) + (EkuU + Ly 


(55) 


where, as we know, U is an arbitrary harmonic solution. Circular orbits are defined by r = r c = constant, so that 
= 0 and, additionally, r, , v | r=ro = 0. From these two conditions it’s possible to evaluate the conserved quantities 

T ,9.-. T 

-r 0 > 0. 

d 2 V I 


|r=ro 0. 

E and L. The orbits will be stable if | r=ro < 0 and unstable if ^^| r=r 

For the sake of saving text space, we do not present here the explicit expression for the derivative (pr |r=r 0 and 
the corresponding values of the quantities E and L. However, in Fig. © we have plotted it for the two first members 
of the family of harmonic solutions Uq and U\ given by (32). The plots show that the second derivative is negative 


for arbitrary ranges of r. So, the orbits of the test particles must be stable in the halo. The plots also reveal that 

the orbits are unstable near the center of the disk. We have also computed these derivatives for other values of the 

angular momentum L , in all cases we have found a similar behavior. Moreover, a simple but long calculation shows 

,2 

that gji- < 0. Thus the particles are attracted towards the center. 


VII. CONCLUDING REMARKS 

We used the formalism presented in (Bj to model an exact relativistic rotating disk surrounded by a magnetized 
halo. The model was obtained by solving the Einstein-Maxwell distributional field equations. In doing so, we 
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introduced an auxiliary harmonic function that determines the functional dependence of the metric components and 
the electromagnetic potential. Accordingly, we separated the total energy-momentum tensor of the system disk-halo. 
Additionally, we expressed the energy momentum tensor of the halo as a sum of two distributional contributions, 
one due to the electromagnetic part and the other due to a material part. As we can see, due to the spacetime here 
considered is non-static (conformastationary), the distributional approach of the Einstein-Maxwell equations allows us 
to work with a strongly non-lineal partial equations system. We considered, for simplicity, the astrophysical consistent 
case in that there is no electric charge on the halo. We obtained that the charge density on the disk is zero. 

In order to analyze the physical content of the energy-momentum tensor of the halo and disk, we projected each 
tensor, in the canonical form, in the comoving frame defined by the local observers tetrad. This analysis allowed us 
to give a complete dynamical description of the system in terms of two parameters (i.e /? and k u ) which determine 
the matter content of the sources. Indeed, the parameter /? in the metric vanishes when it is equal to the isotropic 
pressure and the anisotropic tensor on the material constituting the disk. Similarly, when the parameter k u is equal 
to zero the heat flux on the disk and halo vanishes, a feature of the static systems. So, in this paper we presented, for 
first time, the complete analysis of the most general energy-momentum tensor of a disk-haloes system obtained from 
an exact conformastationary axially symmetric solutions of the Einstein-Maxwell equations. 

The expressions obtained here are the generalization of the obtained for the conformastatic disk-haloes without 
isotropic pressure, stress tensor or heat flow presented in [fjj. Moreover, when we take simultaneously k u = 0 and 
/3 = 1, we obtain its corresponding electrized disk-haloes version. Furthermore, our results are compatibles with the 
description of the relativistic models of perfect fluid disks in a magnetic field presented in HD and the halo presented 
in [Sj. Furthermore, we have shown that the description of the motion of charged particles on the disk and is in 
agreement with the results of analysis of particles motion in the magnetized disks discussed in jllj . In accordance 
with the results presented in Gang, we also shown that the orbits of test particles are stable in the halo for arbitrary 
ranges of r and unstable near the center of the disk. It is also worth noticing that, one can fix the values of the 
parameters /3, k w and b n as well as the number of members of the particular solutions presented here in order to have 
velocity increasing linearly with radius of the disk. 

We have considered specific solutions in which the gravitational and magnetic potential are completely determined 
by a “generalization” of the Kuzmin-disk potential. Accordingly, we have generated relativistic exact solutions for 
magnetized haloes surrounding rotating disks from a Newtonian gravitational potential of a static axisymmetric 
distribution of matter. The solution obtained is asymptotically Minkowskian in general and turns out to be free of 
singularities. 

In short, we concluded that we have presented an exact general relativistic well-behaviored rotating disk surrounded 
by a well-behaviored magnetized halo “material”. In our description we do not impose restriction on the kind 
of “material” constituting the system disk-halo. Consequently, we can speculate that the halo could be made of 
magnetized dark matter. This work provides a solid footing to refine future studies of relativistic disk-haloes systems. 
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Appendix A: The local observers 


We write the metric ([T]) in the form 

ds 2 = -F(dt + ujd<p) 2 + F~ p [dr 2 + dz 2 + r 2 d<p 2 ], (Al) 

where we have rewritten F = e 2 A In addition, we define the tetrad of the local observers //,' y 1 , in which the metric 
has locally the form of Minkowskian metric 


ds 2 = 77 (/i ) (j ,)h {At) <g> h (l/) , (A2) 

h^ t \ = F 1 / 2 { 1,0,0, u;}, (A3a) 

h^ a = F~P/ 2 { 0,1,0,0}, (A3b) 

hW a = F- f3 ' 2 { 0,0,1,0}, (A3c) 

h M a = F~^ 2 { 0, 0 , 0 , r}. (A3d) 

(A3e) 


is given by 
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The dual tetrad reads 


(t) a = A- 1/2 {1 ,0,0,0}, 

(A4a) 

W “ = ^ /2 {0,1,0,0}, 

(A4b) 

•w“ = ^ /2 {0,0,1,0}, 

(A4c) 

F> 3/2 r 

(A4d) 

■(VO - r { 9,0, !}. 


(A4e) 


The circular velocity of the system disk-halo can be modelled by a fluid space-time whose circular velocity V a can 
be written in terms of two Killing vectors t a and ip a , 


v a = v t (t a + jv ), 

where 

u l 11 dt 


(AS) 

(A6) 


is the angular velocity of the fluid as seen by an observer at rest at infinity. The velocity satisfy the normalization 
V a V a = — 1. Accordingly, for the metric (|Al| ) we have 


(vr 


i 

—t a t a — 2f it a ip a — r iip a ip a 1 


(A7) 


with 


t a t a = g tt = ~F 
t Fo. — 9t<p — coF 

<P a <fi a = SW = r 2 F~ p { 1 - A 1+/3 ^), 

consequently we write the velocity as 


FV2 ( i + w n)^i_v? oc ’ 

where 

Vloc = F( + 

is the velocity as measured by the local observers. 


(A8a) 

(A8b) 

(A8c) 

(A9) 

(A10) 
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